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ABSTRACT 



We evaluate the one and two loop contributions to the expectation values of 
^ ' two coincident and gauge invariant scalar bilinears in the theory of massless, 

minimally coupled scalar quantum electrodynamics on a locally de Sitter 
background. One of these bilinears is the product of two covariantly differ- 
entiated scalars, the other is the product of two undifferentiated scalars. The 
computations are done using dimensional regularization and the Schwinger- 
Keldysh formalism. Our results are in perfect agreement with the stochastic 
predictions at this order. 
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1 Introduction 

Quantum field theories which involve either massless, minimally coupled 
(MMC) scalars or gravitons have the propensity for vastly enhanced quan- 
tum effects during inflation. These particles' combination of masslessness 
without classical conformal invariance results in prodigious particle produc- 
tion during inflation [lj . As more and more long wavelength, virtual quanta 
are ripped out of the vacuum, the MMC scalar and graviton field strengths 
grow like the logarithm of the inflationary scale factor. This is evident even 
in free MMC scalar field theory [21 El 13, 



Of course interactions which involve undifferentiated MMC scalars or gravi- 
tons are correspondingly strengthened. 

Powers of infrared logarithms like that in ([1]) arise in the one particle 
irreducible (1PI) functions of a MMC scalar with a quartic self-interaction [51 
El [7]. They occur as well in MMC scalar quantum electrodynamics (SQED) 
El Ell EH] and in massless Yukawa theory [121 02] ■ The 1PI functions of 
pure gravity on de Sitter background show infrared logarithms jTH [151 EE] , as 
do the 1PI functions of Dirac + Einstein [13 [18], and presumably gravity with 
any other theory. Weinberg has recently drawn attention to their appearance 
in fixed-momentum correlation functions [ID, 20J. 

Infrared logarithms are fascinating because they can grow enough during a 
long period of inflation to compensate for even the smallest coupling constant. 
However, what this really means is that perturbation theory breaks down, 
not necessarily that quantum effects become large. Deciding what actually 
happens requires a nonperturbative resummation technique. 

Starobinskii has long argued that the nonperturbative evolution can be 
followed using his stochastic reformulation of inflationary quantum field the- 
ory [21J. Probabilistic representations of inflationary cosmology have been 
much studied in order to understand initial conditions [221 [22] and global 
structure [2H [25] but we wish here to focus on Starobinskii's stochastic for- 
mulation as a wonderfully simple way of recovering secular effects in quantum 
field theory [251 E3 EB] ■ On de Sitter background the technique has recently 
been proven to exactly reproduce the leading infrared logarithms, at each 




(1) 
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order in perturbation theory, for any model of the form [29, 30J, 



(2) 



Provided the potential in ([2]) is bounded below, Starobinskh and Yokoyama 
have used the stochastic technique to give an explicit solution for the late 
time limit |31j . 

Two important generalizations of Starobinskh's technique are necessary: 

• Apply it on de Sitter background to more complicated models which 
also show infrared logarithms, such as Yukawa [12j[T3|, SQED [8| [9| [T0| 
E] and quantum gravity [H EH HQ 13 HI [191 [20] ; and 

• Apply it to a general inflationary and post-inflationary cosmological 
background. 

This first step has already been taken for Yukawa [32] . One of the surprising 
features of the result is that the ultraviolet cannot be ignored, even at leading 
logarithm order. Indeed, leading logarithm results for coincident Green's 
functions can be ultraviolet divergent! This stands in sharp contrast to the 
situation for models of the form (T5]). Aspects of quantum field theory are 
sufficiently counter-intuitive — even in flat space! — that it would be folly 
to ignore the possibility of further surprises as the formalism is generalized 
to models with gauge symmetry and derivative interactions. It is therefore 
imperative to test putative generalizations of Starobinskh's technique against 
explicit perturbative computations at the highest possible loop order. The 
purpose of this paper is to provide this sort of "raw data" for comparison 
with a leading logarithm formulation of SQED [33] . 

We report one and two loop results for the (Bunch-Davies) vacuum ex- 
pectation values (VEV's) of two coincident, gauge invariant bilinears of the 
charged scalar field, 



The Feynman rules are given section 2. Section 3 presents a key result for 
the one loop self-mass-squared, which also serves to fix the divergent parts 
of the conformal and scalar field strength counterterms at order e 2 . The first 
of the two VEV's in ([3]) is evaluated in section 4. Section 5 does the second. 
Our discussion comprises section 6, and the less savory details are consigned 
to an Appendix. 




and 
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2 Feynman Rules 



The purpose of this section is to work out the Feynman rules for SQED 
in de Sitter conformal coordinates. We begin by reviewing the background 
geometry. Then the Lagrangian is given along with a precise definition of 
the renormalization parameters. Of course this allows one to read off the 
interactions in a straightforward manner. We next present the propagators. 
The section closes with a review of the Schwinger-Keldysh formalism which 
adapts the usual Feynman rules to give true expectation values rather than 
in-out amplitudes. 



2.1 Geometry in de Sitter 

We work in the conformal coordinate system of D-dimensional de Sitter 
space, 

ds 2 = a 2 (—di] 2 + dx ■ dx) where a(rj) = — — — . (4) 

Hence the metric is = a 2 r]^ u , where rj^ is the Minkowski metric. The 
affine connection for this background is, 

^ = ffa(fl + ^ + « V ) • ( 5 ) 
For any geometry the scalar d'Alembertian is, 




And the vector d'Alembertian is defined by the relation, 

U;j v = g^f,, pa , (7) 

where a semi-colin stands for covariant differentiation. 

Propagators are most effectively expressed in terms of the length function, 

y(x-x') = aa!{\\x-x* f - {\q-rf\ -i5) 2 } . (8) 

(Here a = a(rj), a' ee a(r/), and the same convention is employed throughout 
the paper. That is, g^ u = g^x) and g' = g^(x').) In the limit that 
5 vanishes, y(x; x') is related to the invariant length £(x; x') between the 
points x M and x ,fl , 

y(x; x) = 4 sin 2 0=Ht{x\ x')) . (9) 
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We often employ the following formula for acting □ upon a function of y(x; x') 
which is analytic everywhere except possibly at y = 0, 



Uf(y) = H 2 {{Ay-y 2 )f"{y) + D(2-y)f(y)} 



V-- 2 f] x ^^-j=5 D {x-x') . (10) 



+Res 



Here Res[F] stands for the residue of F(y); that is, the coefficient of 1/y in 
the Laurent expansion. 

Any scalar which depends upon x M and x ,fl can be considered to be a 
function of y(x;x'). De Sitter invariant vector and tensor functions of x M 
and x ,fl can be represented by including the metric and just the first two 
derivatives of y(x; x') [SJ I3"5] . 



Ha{y8l + 2a'HA Xfl ) , ^ = Ha'(y6° u - 2aHAx u ) , (11) 
J?jL- = H"aa'(y6y u - 2a%HAx v + 2a'HAx,6° u - 2 VfMU ) . (12) 



Here Ax^ = 7]^ v {x — x'Y . Contracting any two of these basis tensors, on 
either primed or unprimed indices, produces a linear combination of the 
basis tensors [3U E5] , 

= H ^- y2) > (13) 

™£e£h = >^~^- < 14 > 

a*>(x)-^ = AH*a' -H 2 -^--^- (15) 

y 1 > dx»dx' v dxPdx'° yva dx ,v dx'°~ 1 ; 

The same is true for covariant differentiation [3U [35] , 

2.2 Renormalizing SQED 

The bare Lagrangian of SQED is, 

£ = -(d^-ie A^ip*(d u + ie A^(pg^y/^ - £ (p*ipRy/^g 

-iAo(pV)V^-JiVWVV=0. (17) 



Here e , £o an d -^o are the bare couplings. No bare mass is required be- 
cause we study massless SQED, and mass is multiplicatively renormalized 
in dimensional regularization. Of course we also study minimally coupled 
SQED, but there is no similar relation for the conformal coupling in dimen- 
sional regularization. Hence £o must appear as a counterterm, even though 
its renormalized value is zero. 

The bare fields are expressed as usual in terms of the renormalized fields, 

if = \fz~2VR and = \JZ 3 A Rlx . (18) 

In terms of the renormalized fields the Lagrangian takes the form, 
C = -Z 2 (d^-i^/z^eoAR^ ip* R (d u + ^fz~ 3 ie A Ru ) tpug* 11 ' y f^-Z 2 ^(pyp R R y f^g 
-^ZlX (^ RVR ) 2 ^g - \z 3 F Rpa F R ^g™^ . (19) 

The various bare coupling constants can be expressed as follows in terms of 
renormalized couplings and renormalization parameters, 

Jz~ 3 e = e + Q , Z 2 f = + 5f and Z 2 2 A = + 5\ . (20) 

Note that we have chosen to make the renormalized 4-scalar coupling zero, 
as we are free to do. Defining the field strength renormalizations as usual, 

Z 2 = 1 + 5Z 2 and Z 3 = 1 + 5Z 3 , (21) 

and dropping the now-redundant subscript R on the renormalized fields, we 
at length reach the form, 

C = -(d vt -ieA l )<p*(d v +ieA v )tpg> a 'J=jj - -F pa F^9 av ^T~9 
-5Z 2 (d lx -ieA lx )(p*(d 1/ +ieA u )(pg^y/^ - 5&*<pRyf=g 

-Ux^^f^- g _ Uz 3 F pa F pv g^g™^r g . (22) 

For the computation reported here we shall not require the interactions pro- 
portional to 5 A and SZ 3 because they do not contribute to either VEV in ([3]) 
at order e 2 . 
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2.3 Propagators 

The scalar propagator obeys, 

y/=g~niA(x; x') = i5 D (x-x') . (23) 

It has long been known that there is no de Sitter invariant solution [36]. If 
one elects to break de Sitter invariance while preserving homogeneity and 
isotropy — this is known as the "E(3)" vacuum [37] — the minimal solution 
takes the form [51E], 



iA(x;x') = A(y) + fcln(aa') where k = g , D , • (24) 

(4tt) 2 T( T ) 

The de Sitter invariant function A(y) is [6], 

A(2/) -(4 7 r)fl f-l V + f-2 V - C0 H 2 J r( f) 



ir(n+D-l)/j/Nn 1 r(n+f + l)/y 



:L n T(n+f) v 4' n-f + 2 r(n- 



•• (25) 



Expression (1251) may seem daunting but it is actually simple to use be- 
cause the infinite sum vanishes in D = 4, and the terms of this sum go like 
higher and higher powers of y(x;x'). Hence the infinite sum can only con- 
tribute when multiplied by a divergent term, and even then only the first 
few terms can contribute. It turns out that most computations in this paper 
require only the derivative, A'(y), expanded to the following order, 

We also need the exact coincidence limit, 

A(0) = -Jfe7TCOt(— ) . (27) 

Hence the coincidence limit of the scalar propagator is, 

lim iA{x; x 1 ) = A(0) + 2k ln(a) . (28) 



6 



The VEV's we seek are both of gauge invariant operators, so it does 
not matter what gauge we use. The photon propagator has been worked 
out in a variety of de Sitter invariant [31] and noninvariant (35] gauges. The 
calculations of this paper happen to simplify greatly in Lorentz gauge, 



d P {V=g<ri 



( a? ;aO} = = ^{ > /-^i 



A, 



(a;; a/)} . 



(29) 



The general form of the photon propagator in any de Sitter invariant gauge 
such as this is 



or \ 9 2 y dy dy 

B (y) + c (y> 



dx^dx'' 



dx^ dx'' 



(30) 



For Lorentz gauge (1291) the two functions B(y) and C(y) can be expressed as 
follows in terms of a single function 7(2/) [34"), 138] . 



B(y) 
C(y) 



1 



A(D-1)H 2 
1 



A(D-1)H 2 



(4y-y 2 h'(y) - (D-l)(2-yh(y) 



+{2-y)i{y)-{D-l) 1 (y) 



(31) 
(32) 



The gauge condition fl29|) is obeyed for any function 7(1/). What fixes 
7(7/) is the equation for the photon propagator [38] . 



9W-R U 



(x; x') = ig flv S D (x—x) + v ^gd^d' L ,iA(x; x') . (33) 



The unique de Sitter invariant solution for 7(7/) is J3S] 



* (4vr) 2 I "2 



r(f)/4^-i 



any 

■1 V 



+ 



D-3 



E 

n=0 L 



r(n+D-l) /yxn 



r(n+f + 



x(n+l)[v>(2-— -l) +^(n+D-l)-^(n+2 



r(n+f +1) /ysn+a-f 



(!) 



T(n + 3) M 

x 



n+3- 



2 



0(2-^)-V«(y-l) +V(n+y + l)-V(n+4-y) 



Here the symbol u ip(z) n stands for the polygamma function, 

^) = ^ln(r(z)) . 



(34) 



(35) 



As with the scalar propagator, we do not require the full complexity of 
7(7/). For most computations in this paper the following expansions suffice, 



i{y) 



i{y) 



B{y) 



C{y) 



1, s ,Ds H D - 2 { 2 /4Nf-i 



HI 



+ — + 0(D-4) , (36) 



± 4 



-(D-1)T 



D^ H°- 2 f/4 N £ 



(47r)^" 
21n( 



£-2— ?- 
2 L>-2 



:i-D 3 



4 V 2^( 47 r)# 1 £>-2Vy 
ln(|) 



16 v 2 /( 4 ^)f 

2 Ml) 



(l-I) 2 

.(!)#_ 

\ y j 
21n(|) 



:i-D 2 



2 MS) 



ft*" 1 

0(D-4) 



1-1 



+ 



4n#-i 



(l-I) 3 (l-I) 2 ( 
We also need the exact coincidence limit, 
^ 0) = l ( D-l,H^T(D-l) 



D-2 



2VD-3/ (47r)f T 



X 



(37) 



0(£>-4) , (38) 



(39) 



y(2-^)-V»(y-l) +^-l)-V(2)]|. (40) 



The coincidence limit of the photon propagator follows from this, combined 
with expressions flU]), (J3U132D and (PIE2]), 



lim z 



1 OC - <3j J 



7(0)5 



(41) 



2.4 The Schwinger-Keldysh formalism 

The Feynman rules we have just presented would suffice for computing the 
matrix element of any operator between the state which is free vacuum in 
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the asymptotic past and the state which is free vacuum in the asymptotic 
future. These in-out matrix elements are well adapted to scattering experi- 
ments in flat space but they do not correspond to observations that can be 
performed in de Sitter. There is no S-matrix in de Sitter [39, 40J. In fact the 
vast expansion of spacetime in the infinite future means that in-out matrix 
elements even diverge off-shell [41J. 

The physics reason behind the math problem is that free vacuum is an 
infinitely poor guess for the state after an endless history of inflationary 
particle production. Indeed, the fact that we do not know what becomes 
of the state in the infinite future is the whole reason the computation is 
interesting! Under these circumstances, a more realistic theoretical question 
is, what happens if universe is released in free vacuum (or some other state) 
at a fixed time and then left to evolve as it will? One can follow this evolution 
by computing the expectation values of time-dependent, local operators in 
the presence of the Heisenberg state. 

Of course expectation values of the sort just described could be com- 
puted canonically. However, the canonical formalism is extremely cumber- 
some because it does not reflect the underlying spacetime symmetries in a 
manifest way. Schwinger long ago modified the ordinary Feynman rules to 
provide a manifestly covariant procedure for computing expectation values 
[II, S3 IH EH] • Many excellent reviews of this formalism exist gSj 1471 1381 [T9] 
so we will merely give the key identity which relates it to the operator for- 
malism 



r*(<%M)r(0!M)|*) 



[#+][#_] 5 <j>^{£)-<j> + {l) 



[0_]O 1 [0 + ]^[0_( S )]e l/ " d '{ L[ * +W1 - L[0 " W1 }vl>[0 + ( s )] . (42) 



Expression fj4*2l) is formulated in the context of a scalar field <p(x) whose 
Lagrangian (the space integral of the Lagrange density) at time t is L[(f)(t)]. 
The left hand side of (I42p is the canonical expectation value of the product 
of two arbitrary operators, 0\ [(f)] and O2 [(f)] , in the presence of a Heisenberg 
state |\&) whose wave functional in terms of the fields at the starting time t = 
s is \l/[0(s)]. We have time-ordered Oi[(f>] and anti-time-ordered C?2 [</>]■ The 
right hand side gives the functional integral representation of this canonical 
expectation value. The last time t — i can be any point in the future of the 
latest operator in either Oi[(f>] or 02 [</>]■ 
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It is important to remember that there is only one type of operator for 
each of the fundamental fields. The (D-number functions (f>± which appear 
on the right hand side of fj42|) are merely dummy variables in a convenient 
functional representation for the canonical expectation value on the left hand 
side. From this functional representation we can read off the following rules: 

• Each endpoint of each propagator carries a "polarity" which can be 
either "+" or 

• External lines emanating from operators in Oi[<f>] are "+", whereas 
those from operators in 2 [0] are "— " ; 

• Any given vertex is either all "+" or all "— " ; 

• The "+" vertices are identical to those of the in-out formalism, whereas 
the "— " vertices have the opposite sign; 

• The state wave functional \l/[0(s)] is assumed to be free vacuum, plus 
possible perturbative corrections which would show up as interactions 
on the initial value surface; and 

• For SQED in de Sitter, the various polarities of the propagators are 
obtained from the Feynman propagators of the previous subsection by 
making the following replacements for the length function y(x; x') [oD] , 



It is worth calling attention to some consequences of these rules, 

• To each iV-point 1PI function of the in-out formalism there correspond 
2 N 1PI N-point functions in the Schwinger-Keldysh formalism; 

• The absence of mixed polarity counterterms means that no mixed po- 
larity 1PI functions harbor primitive divergences |49j; and 

• The propagators of mixed polarity obey the homogeneous analogues of 




y++(x; x') 
y+-(x;x') 




(43) 
(44) 
(45) 
(46) 



(J23D and ([33]). 
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3 One Loop Self-Mass-Squared 



The purpose of this section is to derive a key result for the one loop self- 
mass-squared. By renormalizing this quantity we also determine the order 
e 2 contributions to the renormalization parameters 5£ and 5Z2- 

The three diagrams which contribute to the ++ polarization of the one 
loop self-mass-squared are depicted in Fig. 1. 




Fig. 1: One loop contributions to —iM 2 (x;x'). 

Of course the left-hand diagram is the most difficult. In an arbitrary gauge 
it would consist of four terms. However, one of the nice things about Lorentz 
gauge ( |29i) is that each of these terms makes the same contribution, 



v-gg 



I Ivu ■ 

-9 9 « 



-gg 



hp 



-g'g wa * 



-gg 



tip 



A a (x;x')d IJ ,d > u iA(x;x') 

(x; x')& v iA(x\ x')\ 
(x; x')d fl iA(x; 
x; x')iA(x; 



-9'9 ,va i 



-g'g' va i 



P A a 



= -Ae 2 ^g^^g~'g'™i 
The next step is to take the derivatives, 



{x;x')d„dliA(x;x') . (47) 



d„&JA(x;x>) = A'(y)^ + A"{y)^^ + 5^-^5 D (x-x') , (48) 



Q x HQh 



dx p dx 1. 



and carry out the tensor contractions with the aid of identities (I13H15I) . One 
then substitutes relations (13T1) and (1321) . The analysis is, 



-AeW-gg^xl-g'g'^i 



P A a 



[x; x')d fl d' I/ iA(x; x) 



d 2 y 



dy dy 



-4e 2 V=gg" p \/-g'g' ua \ + c^-^ faffjAfa x') , (49) 



dx p dx' a 



dx p dx' a 
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-Ae 2 H 4 ^J-g'l(Ay-y 2 )A" (2-y)B + (4y-y 2 )C + ADBA' 



(4y-y 2 )A' B - (2-y)C \ + iAe 2 y{0)^5 D (x-x') , (50) 



= Ae 2 H 2 ^^-g'^y-y 2 )[A" 1 + A'i\ + £>(2-y)A' 7 j 

+i4e 2 y(0) v ^g~5 D (x-x) . (51) 

By comparing with (|T0|) we can recognize (151!) as the scalar d'Alembertian 
acting upon the indefinite integral of A' (y)y(y). Because we shall often need 
to take the indefinite integral of different functions F(y), we shall denote it 
with the symbol I[F] as follows, 



I[F](y) = J dy'F(y') . 
By consulting the expansions (}2"E]) and (jMj) we see that, 

ttD-2 



Res 



H D ~ 2 „sD 



4ttt ^ 2 



r --i h(o) 



(52) 



(53) 



Therefore the final term in (I5TT) gives precisely right residue contribution and 
we conclude that the first diagram of Fig. 1 can be written as, 



-4e* 



'99 



hp 



-9'9' v °i 



(x; x')d,d' u iA(x; x') =Ae 2 ^J -g'UI[A'i\ . 



(54) 



The other diagrams of Fig. 1 are straightforward. The middle one is, 



-ie 2 ^g\g^i 



(x;x')8 D (x-x f ) = -ie 2 Dy{0)y/=g5 D {x-x') . (55) 



The rightmost diagram of Fig. 1 represents field strength and conformal 
renormalizations, 



i6Z 2 y/=gn6 D {x-x') - i5£R^/=g~8 D {x-x') . 



(56) 



Adding f|54l) . (|55|) and (1561) gives our result for the ++ polarization of the 
regulated one loop self-mass-squared, 



iM 2 ++ (x;x r ) = ^^gJ-g'U^A'y] + i5Z 2y ^gU5 D (x-x' 



—i 



D(D-l)H 2 8£ + e 2 Dy(0) 



(x-x') . (57) 
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To renormalize the one loop self-mass-squared we must isolate the diver- 
gences implicit in the first term of (157)) . From the expansions (125]) for A'(y) 
and (136)) for 7(7/) we compute, 



A\vWv) = -kD-i)r(j)r(j-i 



H 



2D-4 



(4tt) 



/4\£>-i /4\2 /4 
- +3(- +3(- 

\yj \yj \y 



+ 



/4n2 ,4x 4 3 



1-1 



D-4 

2/ 2 



(58) 



Hence the indefinite integral is, 




It is important to understand that one typically wants to use —iM 2 



[x: x 



inside an integral over so an expression is "renormalized" when it has 
been written in a form which is integrable (with respect to x ffl ) in D = 4 
spacetime dimensions. Each term of (159)) except the first meets this require- 
ment. To isolate the divergence in this first term we exploit relation (TTUT) to 
express it in terms of a lower power which integrable, 



(T 2 



y 



(D-3)(D-A) H 2 




D-A^y 



/4\-D-3 



(60) 



The d'Alembertian with respect to x M can be pulled outside the integra- 
tion over x ,fl , so (160)) really is an integrable expression. We could take the 
limit D — > 4 at this point, were it not for the factors of 1/(D— 4). Of course 
these represent the ultraviolet divergence we have been laboring to extract! 
To isolate it on a local term suitable for renormalization, we add zero in the 
form, 



¥■ TT-D 







(AirpH 



5 D {x-x') 



r(f 



1 




Hence we conclude, 
/4nD-2 2 

\y) 



f(47 



n)fH- D id D (x-x') 



(61) 



(D-3)(D-A)\ r(f-l 



H 2 



V 



\yj 
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-AO 



4xD- 3 D{D-2) /4xf-i 



D-A\y y 



ZlT-D 



(Att)-H 



S D (x-x r ) 



(D-3)(D-A) r(f-l) 



8(£>-3) ^ 



, (62) 



H 2 [y M 



+2 d)K!)-G)+°( fl - 4 >- «b) 



■y 



Substituting fl63l) in expression fl59l) . and the result into fl57|) . gives a 
form for the self-mass-squared from which renormalization parameters can 
be inferred, 



-iMl + {x;x') = ihz 2 + 



e 2 H D - 4 2(D-l)T{%-l) 
(4tt) 4 (D-3)(£>-4) 



—i 



D{D-l)H 2 5£ + e z D~t(ti) ^)5 D (x-x') - 



3e 2 H 2 
(Atx) 1 



x 



° 2 H!) + 



3e 2 H A 
(4tt) 



1 /y 



n=l 



+ 



4( ! ) _ 41n (|) + 8K 1 -|)- I 4 ? ]ln(|)} + Op-4).(64) 



This fixes the one loop counterterms up to finite renormalizations, 

e 2 7 (0) ^ 12e 2 ^ fin 



SZ 2 



+ 



(D-1)H 2 D(D-l) 



+ 0(e 4 ) 



(65) 



e 2 ^" 4 2(D-l)r(f-l) 

+ e dZfin + 0(e ) 



( 47 r)f (£>-3)(D-4) 
Taking the limit .D — > 4 gives the renormalized result at one loop order. 
-iM 2 ++ (x; x') = ie 2 5Z^^^U5\x-x') - H2e 2 H 2 5^ n ^g5 4 (x-x') 



3e 2 H 2 
' (4tt) 



3e 2 ff 4 
(4tt) 



71=1 " ^ 



4 _ _4i n * +8 i„ h_£ 



6 



ln(|j • (67) 
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4 (Q\(D^)*D u tp\tt) 

The purpose of this section is to evaluate the leftmost of the VEV's in ([3]) 
at one and two loop orders. We begin by decomposing the gauge invariant 
operator product into noninvariant components, 

{D^)*D u y = (d^-ieA^ip*(d v +ieA v )ip , (68) 
= d^ip*d v (f - ie(f*A^d u (f + ied^*A u Lp + e 2 ip*A fl A u tp . (69) 

The next three sub-sections are devoted, respectively, to writing out and 
partially evaluating the expectation values of the order e°, e 1 and e 2 operators 
in (169|) . These results are combined in the final sub-section, at which point an 
important cancellation occurs. We then give explicit results for the remaining 
integrals. 

4.i {n\d^*dMty 

The only diagram which contributes at one loop is Fig. 2. 




Fig. 2: One loop contribution to (ft\d fl <f*(x)d v ip(x)\ft) . 

It corresponds to the coincidence limit of the differentiated propagator, whose 
value comes entirely from the n = 1 term in (f25|) . 

lim d.d'jAix; x') = H ° I j^ff - x --H 2 9liu = -(^^)kH 2 g. u . (70) 
x'^x v ' ; ( 47r )# r(f + l) 2 y ^ V d J ui»> \ ) 

Fig. 3 depicts the three diagrams which contribute at two loop order, the 
scalar bilinear being an insertion at point x^. 
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X 

a" 



X X 



x'«WVWW^x" -j- T i + T i 



Fig. 3: Two loop contributions to (Q\d^ip*(x)d u ip(x)\Q) . 

The alert reader will recognize these as the three diagrams of Fig. 2, with 
the external lines joined. This means we can express the diagrams in Fig. 3 
as a sum over simple integrals involving the four one loop contributions to 
the self-mass-squared, 

]T f d D x'd^iA +± (x;x') f d D x"djA +± (x;x") x -iM 2 ±± (x';x") . (71) 
±± J J 

For our purposes it is better not to use the renormalized self-mass-squared 
(j67j) . which is only valid inside integrals over suitably smooth functions and 
in the limit D — > 4. It is instead superior to use the exact, regulated result 
( I57p . Symmetrizing the d'Alembertian and taking account of our result (1651) 
for the conformal counterterm, we can write the ++ polarity as, 



-iM 2 ++ {x';x") = 2e 2 ^ s /^(n' + n")I[A' 1 }(y, 



+iSZ 2 ^/-g'n'5 D {x'-x") -H2e 2 H 2 S£ &n ^-g'5 D {x-x') . (72) 
The other polarizations give, 

-iM 2 + _(x';x") = -2e 2 f^>f^< (□' + n")/[A' 7 ](y + _(x';x")) , (73) 
-tM 2 _ + (x';x") = -2e 2 ^^{n' + U")I[A'l}(y- + {x'^")) , (74) 
-iM 2 __(x';x") = 2e 2 ^^'{U' + n")I[A' 1 ](y_ 



-i5Z 2 ^g~ f D'S D (x'-x")+il2e 2 H 2 6^ &n ^g~ f 6 D (x-x') . (75) 

The counterterms are simple on account of the delta functions. The 
contribution to (|71|) from the conformal counterterm is, 

-il2e 2 H 2 5^ n J d D x'^g~'^iA ++ (x; x')djA ++ (x; x ) 

-0„zA + _ (x; x')d v iA + _ (x; x') ) . (76) 
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We postpone evaluation of these integral to the end of this section. On the 
other hand, the contribution to (ITTj) from the field strength counterterm can 
be evaluated right away, 




— d fl iA + _(x; x')d u n'iA + _(x; x') 
= -5Z 2 lim d^&JA^x; x') = 5Z 2 (^±)kH 2 g^ . (77) 

x'^x v J_J J 

The nonlocal contributions suggest a partial integration, which raises the 
issue of surface terms. In the Schwinger-Keldysh formalism the various po- 
larizations combine to completely cancel contributions from the future time 
surface, and from the spatial surfaces [50J. However, there are nonzero con- 
tributions from the initial value surface. We shall assume that order e 2 
corrections to the state wave functional completely cancel these surface con- 
tributions. Even if this is not correct, the surface contributions should fall 
off exponentially [6]. 

Of course the point of partially integrating is to take advantage of the 
simple results of acting d'Alembertians on the various propagators, 

sf^g'U'iA ++ {x] x') = iS D (x-x') , sf^g'U'iA + _{x] x') = , (78) 

y^D"2A ++ (x;x") = t6 D (x-x") , yf^g 7 'D"iA + _(x; x") = . (79) 

One should also note that the derivative acting upon whichever outer leg 
propagator contributes a delta function, will ultimately wind up acting to 
undo the indefinite integral of A'j, 

Jd D x"^d u n"tA ++ (x; z")/[A' 7 ](Wa/; x")) =d l/ I[A' 1 ](y ±+ (x'; x)) , (80) 

= ^-A'(y +± (x;x')) 1 (y +± (x;x')) . (81) 

Hence the nonlocal terms give, 

z4e 2 / ^x'^^d^A^x^A'iy.My^) - (+-)} . (82) 

(Parenthesized indices are symmetrized throughout this paper.) We shall 
not evaluate this contribution because it cancels against one in the next 
sub- sect ion. 
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4.2 -ie(n\(p*(x)A fi (x)d„(p(x)\n) + ie{Q,\d^*(x)A v (x)(p(x)\Sl) 
Both VEV's have the diagram topology shown in Fig. 3. 



x 



x 

Fig. 4: Lowest contribution to 

-ie(n\ip*(x)A II (x)d u ip(x)\n) + ie(n\d fl ip*(x)A u (x)ip(x)\Q). 

The leftmost VEV can be simplified by another partial integration whose 
surface term on the initial value surface we shall again assume is canceled by 
an order e 2 correction to the state wave functional, 



le 



2 d D x'J-g'g^li 



Ap 



& a iA ++ (x;x')djA 
—iA ++ (x; x')d u d' a iA ++ (x; x) 



+- 



i2e 2 d D x'J-g'g 



Ap 



fx; x')d' a iA ++ (x- x')d u iA ++ (x; x') - (+-) . (83) 



The rightmost VEV makes the same contribution with the indices /i and v 
interchanged. 

The next step is to act the d' a , carry out the contraction using identities 
( TI3HT4"]) . and substitute (I31H32I) . Because none of these operations depends 
upon the Schwinger-Keldysh polarities we will carry them out generically, 



g'o°i 



Ap 



[x; x')& a iA(x; x') 



dx^dx'f 



dx^ dx' p 



dx h 



3 V Al I 



dxt 1 



A\y) (2~y)B(y) + (Ay-y 2 )C(y) 



kH . 






Ao 



dx^ 



A'(y)j(y) 



kH r 



Ao 



x: x 



(84) 

(85) 
(86) 
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Therefore the two VEV's of this sub-section contribute, 



-iAe 2 j d D x'^-^{d u) iA ++ (x-,x')A'(y ++ My ++ ) - (+-)} 

-i4e 2 kH J d D x'a' D - 1 {d {ll tA ++ (x;x , )t[ u) A \(x;x') - (+-)} . (87) 



Of course the first line of (1571) just cancels (J82J), so we need not evaluate 
either of them. We postpone evaluating the second term until the end of this 
section. 



4.3 e 2 (n\(p*A ll A l/ (p\n) 

The only order e 2 diagram that contributes to this VEV is depicted in Fig. 5. 



Fig. 5: Lowest contribution to e 2 (fl\ip*(x)A lJi (x)A u (x)ip(x)\Q) . 
It is a straightforward coincidence limit, 

(x;x)iA(x;x) = e 2 7 (0){A(0) + 2Hn(a)}^ . 



4A The Final Result 

Adding the various contributions — (}70]h (1761) . (1771) . (152|). (1571) and (fs 
gives the following result for the first VEV, 



n 



■D-l 



(D^)W uV \n) = -[l-5Z 2 ](—)kH 2 g^ + e 2 1 (0){A(0) + 2k\n(a)}g lliy 



Q2e 2 /7 2 ^ fi n/^a;' v /^ ; {^A ++ (x; x')djA ++ (x; x') -(+-)} 



-iAe 2 kHjd D x'a' D ' 1 {d ifl tA, 



(x:x )i 



v)Aq 



(x;x')-(+-)} + 0(e 4 ).(89) 
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It remains to evaluate the two final terms of (1891) . For that purpose the 
2nd covariant derivative of a function f(y), 



UM = H 2 (2-y)f\y)g, v + f"(y) 
implies an identity we shall employ many times, 



dy dy 
' dx^ dx v 



(90) 



F(y) 



dy dy 
dx^ dx v 



-H 2 g, u {2-y)I[F}{y) + d,d u -T? u d p I 2 [F}(y) . (91) 



In using ( 19TT) . one extracts derivatives with respect to x M from integrations 
over For example, consider the integral associated with the conformal 
counterterm in (189~|) . 



d D x'^{d ll tA ++ (x; x')d u iA ++ {x; x') - (+-) } 

d D x'yPg~'L' (y ++ )^ + kHa6°] [A {y ++ ) $L + k HaS° u 



- +- 



(92) 
(93) 



-H 2 g, v Jd u x^-g> (2-y ++ )I[A' 2 ](y ++ ) - (+-) 
+ 



d^-T^d, \d u x'^-g' I 2 [A'%y 



+- 



+2kHa^d v) d u x^-g> A(y ++ ) - (+-) , (94) 



= -H 2 gflu Jd^x'^-g'[(2-y ++ )I[A' 2 }(y ++ ) - (+-) 
+H 2 {-g, u + a 2 5l5l[d laa ~l]}d Xaa fd D x'^'[l 2 [A 2 ]{y 



+2kH 2 a 2 5Xd lna ld n x'J-g> A{y ++ ) - ( 



D i 



(95) 



Here d\ na stands for the derivative with respect to ln(a), which is the only 
coordinate upon which the integrals can depend. 

Reducing the final term of (I59"j) to a similar form is facilitated by two 
relations for temporal basis tensors, 



^L = Ha'(y-2) + 2Ha and — 



Ha'p- + 2H 2 a 2 5°„. (96) 
dx u v \ i 
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Combining these with (I3"0"j) and (I31ti32j) . we obtain the relation, 



A (x;x) = B(y) 



d 2 y 



C 



dy dy 



Ha 



dx v drj' dx v drj' 



(97) 



dy 



B{y)-{2-y)C{y)\ + 2Ha-^C(y) + 2H 2 a 2 ^B{y) , (98) 



dx u 

Substituting this relation into the final term of ( 1591 . decomposing the deriva- 
tive of the scalar propagator, and making use of (19T|) results in six integrals 
whose sum yields the final term of ( |89|) . 



<9y <9y 



DllD 



'-{A'(y ++ W(y ++ ) -(+-)} 



(100) 



--^^|^y^a;V-[(2-y ++ )/[AV](27 ++ ) - (+-) 



+ 



iSe^kH'a fd D xW D -^^{A'( y++ )C(y, 



j, (101) 
-(+-)}, (102) 



i8e 2 H A ka\g^ld D x'a' D - 1 (2-y ++ )I[A'C](y 



+ 



J 2 L4'C](y ++ ) - 



5f^-a 2 5j5°(9i na -l) ^na/d 1 

J^(lna) = -z8e 2 ^ 3 a 2 /dVa^- 1 5^{A'(y ++ )5(y ++ ) - (+- 

= -i8e 2 ifefr 4 oa 2 <5j^a tao ydVa /D - 1 {/[A / S](y ++ ) - (+-)}, 
/^(ln a) 



iAe 2 k 2 Ha 



D-l 

iAe 2 k 2 H 2 



^fa^Xd^J^x'a^^) - (+-)} , 



/J,(lna) ^ -i8eWa 2 fd D x>a' D -%^{C(y 



- +- 



= -?8e 2 A; 2 J ff 4 aa 2 (5j^^ na yrf D x , a /D - 1 {/[C , ](i/ ++ ) - (+-)} 
/^(lna) = -z8e 2 A; 2 # 4 aa 2 c^° Jd D x'a ,D - 1 {B(y ++ ) - (+-)} 



,(103) 

)},(104) 
(105) 

(106) 

(107) 

(108) 
(109) 
(110) 
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The next step is compute the various functionals of propagator functions 
which appear in (j95p and Ij^ 6 - Of course we take D = 4 on any terms which 
make finite contributions. For example, the square of A'(y) is, 

Its integral is, 



4 {4tt) d { D-lV D-2V K V 

(112) 

Hence the integrands of the first two terms in (}95i) are, 



+16 + 0(^-4)1, (113) 



2 



-4\ D-2 



+2(i)-41n(|)+0(C-4)|. (114) 



The integrand for the final term in fl95|) is just the D = 4 limit of A(y), 

H 2 



Mv) 



(4^{$- 21n (f)- 1 + ^- 4 )}- ( 115 ) 



It turns out that each of the 11 integrands required for evaluating ( l95l) 
and J*" 6 can be expanded in terms of just two potentially divergent functions 
of y/4 and eight finite ones. The two potentially divergent functions are, 

A\D-1 j /4\D-2 

(y) aM (y) ■ ^ 

In considering the finite functions it conserves space to set x = y/4. These 
functions come from multiplying the less singular parts of propagators, de- 
composing by partial fractions, and then performing whatever integrations 
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x(l- 


-x) 


1 




x(l- 


x) 2 


1 




x(l- 


x) 3 


1 




x 2 (l- 


-s) 


1 




x 2 (l- 


x) 2 


1 




x 2 (l- 


x) 3 



and subsequent multiplications are prescribed. The partial fractions decom- 
positions are facilitated with the familiar identities, 

1 1 1 

= - + T— ' 117 

x 1—x 

1 1 1 

x 1—x (1 — x) 2 ' 

= - 2 +- + 7^-> (120) 
x z a; 1 — x 

12 2 1 . . 

x z x 1—x (1— x) 2 
1 3 3 2 1 

The integrals we require are, 

/dx^ = (123) 

J X 2 XX 

Jdx^l = \ln\x), (124) 

J dxln(x) = xln(x)— x, (125) 

J dx\n 2 (x) = xln 2 (x)-2xln(x) + 2x , (126) 

dx , V ; Q = : V \ — In x ) + -ln 1-x - — -, 127 

(1-x) 3 2(l-x) 2 2 v ; 2 v ; 2(l-a;) V ; 

= ^-ln(x)+ln(l-x), (128) 
(1 — x) 2 1 — x 

[ d J^L = -l^i-^ln^-g^, (129) 
i 1-x ^n 2 

J dx\n(l — x)\n(x) = — (1— x) ln(l — x) ln(x) + (l — x) ln(x) 

+ y 777 • (130) 

n t>(n + l) 2 1 ' 

Contributions which are analytic at y = — such as the infinite sums - 
cancel when we take the difference of ++ and H — terms. So the eight finite 
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Integrand 



Prefactor 



'4\D-1 

1 y 



•4\D-2 



(2-y)I[A' 2 ](y) 



F 2 (f)H 2 



4(477)° 



D-l 



2 + 4 



L>-2 1 D-l 



I 2 [A' 2 ](y) 



F 2 (f)H 2 



(4tt) d 



1 



(D-l)(D-2) 



(2-y)I[A'Y](y) 



(D-1)V 2 (^) H 2 
8(4tt) 1 



D-l 



-2 + T^T + ' 



D-l 1 (D-2) 2 



/ 2 [A' 7 '](y) 



(D-i)r 2 (f ) g 2 

2(4tt) 1 



1 



(D-l)(D-2) 



(2-y)I[A'C](y) 



Y 2 (%)H 



16(4?r) 



2D-6 
ZS 



D-l 



2D 



D-l D-2 (D-2) 2 



I 2 [A'C](y) 



T 2 (f)H 2 



4(4tt) d 



1 



(D-l)(D-2) 



I[A'B](y) 



T 2 (f )H 2 



4(4tt) d 



(D-2) 2 



Table 1: Potentially Divergent Contributions to Expansions of Integrands. 



functions can be taken to be, 



fi(x) 



fs(x) = 



X 

ln(x) 

x 

ln(x) 
(1-x) 2 
ln(x) 
l — x 



f 5 (x) =ln 2 (x) -21n(l-x)ln(x) , 
fe(x) = xln 2 (x) — 2xln(l— x) \n(x) 
fr(x) = ln(ar) , 
/ 8 (x)=ln 2 (x). 



(131) 

21n(x), (132) 
(133) 
(134) 



We have tabulated the expansions for the integrands in (}95l) . and those for in- 
tegrals /*~ 6 , in Tables [T] and [2J Tabled] gives the numerical coefficients of the 
potentially divergent contributions and Table Ogives the finite contributions. 

The next step is to perform the integrations, against either a' D or a /D_1 , as 
the case may be. We have done this generically for each of the two potentially 
divergent functions in Table [3l Table H] gives the integrals of the eight finite 
functions times a' 4 . The analogous results for a' 3 are reported in Table EJ 
The Appendix describes how these results were obtained. 

One evaluates each of the 11 integrals in (|95|) and J^~ 6 by multiplying 
the appropriate tabulated results, and then acting whatever derivatives are 
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Integrand 


Prefactor 


fl 


h 


h 


h 


h 


h 


fi 


(2~y)I[A' 2 ](y) 


H 4 

A ( Att\4 
4^47T ) 























I 2 [A' 2 ](y) 


H 4 


2 

















-4 


My) 


H 2 


1 

















-2 


{2-y)I[A'i](y) 


3H 4 
o^47r j 


13 


4 


6 


4 


-10 


20 





I 2 [A'y](y) 


3H 4 
2(4tt) 4 


i 

2 








-3 


1 


-5 





(2-y)I[A'C](y) 


H 2 
16(4tt) 4 


-20 


-8 


-6 


-16 


18 


-36 





I 2 [A'C](y) 


H 2 
4(4vr) 4 


5 
2 








3 


-2 


9 





I[A'B](y) 


H 2 
4(4vr) 4 


-6 


-3 





3 


5 








i{y) 


3H 2 
2(4tt) 2 


1 





1 














i[C](y) 


1 

4(4tt) 2 


1 





-1 


-2 











B(y) 


1 

4(4tt) 2 


-1 





-1 


-2 












Table 2: Finite Contributions to Expansions of Integrands. 



/(*) 


Jd D x'a' D {f(^)-f(^)} 






(D-2)(D-4) +°(„' D - 4 ) 


+ O(^,D- 4 ) 


e°- 2 


7^4+00-4) 


(D-3)(D-4) + (a2' D - 4 ) 



Table 3: Divergent Integrals. Multiply each term by 



25 



/(*) 




1 


-5 + °(s) 


ln(x) 

X 




ln(x) 
(l-x) 2 


+ T + 


ln(x) 
l-x 




ln 2 (x)-21n(l-x)ln(x) 


1 vr 2 , /l/ln(o) \ 
6 9 V a > 


x ln 2 (x)-2x ln(l-x) ln(x)-2 ln(x) 


5 vr 2 i >o/ln(o)\ 
24 18 "'" a ) 


ln(x) 


|ln(a)-| + 0(i) 


ln 2 (x) 


|l„ 2 (a) - |ln(.) + \ ~ f + O(^) 



Table 4: Finite de Sitter Invariant Integrals. Multiply each term by 





/dVa«{/(M)-/(t)} 


i 

X 


-ln(a) + 1 + 0{\) 


ln(x) 

X 


±ln 2 (a) + 21„(a) 3 + f + O ( ^ ) 


ln(x) 
(l-x) 2 


ln(«) - 1 - £ + O(J) 


ln(x) 
l-x 


iln 2 (a)-21n(a) + 3+0(^) 


ln 2 (x)-21n(l-x) ln(x) 


— ln 2 (a) + 21n(a) — 5 + £ + O(^) 


x In 2 (x)-2x ln(l-x) ln(x)-2 ln(x) 


_I ]na(a) + | ln(a) _g + ^ + (3^1) 


ln(x) 


I a -ii n(a ) + i + 0(I) 


ln 2 (x) 


|aln(a)-^-a-|ln 2 (a) + 21n(a)-| + £ + 0(^) 



Table 5: Finite de Sitter Breaking Integrals. Multiply each term by 
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prescribed. For example, the first term in (1951) gives, 

f 2 



-H 2 g^ x \ 2 ' - _ x 



4 (4tt; 



+ 



# D r(f-l)\ D-l (D-2)(D-A) 



4 



D-2 



r(f; 



D-4 



( 47 r)f 2(D-4) 



D-2 L>-1 

The second term of (|95|) gives an even simpler result, 

D \ H 2D ~ A *(4vr 



fiV • (135) 



H 2 {- giw + a 2 5l5 Q u [d lna -l]}d Xna x r 
f 1 



x 



2^ (4vr) D " H D T(%-1) 



X' 



2 1 

— x h 2 x 4 x 

(D-l)(£>-2) D-4 2 



- ma 

6 K J 36 



iH 2 



247T 2 



The final term of f l95l) is finite from the start 

1 



F 2 ?'lfi7T 2 



167T 2 



2x 



-ma 

6 v ; 36 



(136) 



iH 2 

—— a 2 5°5° 



(137) 

Summing the three terms and multiplying by —H2e 2 H 2 5C,fi n gives the total 
contribution from the conformal counterterm, 



(4tt; 



6 



D-A 



(138) 



In these and all subsequent formulae we neglect terms which vanish at D = 4, 
or which fall off like 1/a. 

The procedure for J^~ 6 is the same. Consulting the appropriate tables, 
acting derivatives and summing, we find, 



I 1 = 


e 2 R 2D-A 




I 2 = 

flU 


e 2 H 2D-4 






+ 


-(D- 



(-D 2 + 6D-6)T(D-V 
(D-2)(D-A) 



-9 + vr 2 Sg^, 



(139) 



(£> 3 -5-D 2 + 14D-14)r(.D-l) 2 

f7T 



(D-1)(D-2)(D-3)(D-A) 



(D-1)(D-3)(D-A) 



81n 2 (a) + 321n(a)-58 + -vr 2 aW! ,(140) 
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{4tt) d { (£>-2)(£>-3)p-4) 



+81n 2 (a)-321n(a) + 68-V aW , (141) 



r =o 



2 



J 5 



e 2 H 



(47TJ 



e 2 R 2D-A 



8 In' (a) - 32 ln(a) +48 - -tt 2 j a 2 ^, 
-8 ln 2 (a) +32 ln(a) -48 + |tt 2 | a 2 ^<5°. 



(4tt) d 

The sum of all six integrals is, 

-D 3 + 9D 2 -20L> + 16)r(D-i; 



(142) 
(143) 

(144) 



6 2o-2D-4 

Ejk ( " 

fc=l 



(47T) 



(D-l)(£>-3)(D-4) 

A(2D-3)T(D-1) 



-9 + 2tT 



+ 



(L>-l)(L>-2)(D-3)(L>-4) 
Note the cancellation of all infrared logarithms. 



+ 10 aW .(145) 



We can now substitute 



n 



1-6 Zn 



and (I145P in (18"9"|) to give our final result, 
D-l- 



D 



^V + e 2 7 (0){A(0) + 2A;ln(a)}^ 



fin 



(4tT 



6 



D-4 



^ - 8 a W + 



e 2 H 2D-A 



(4+ 



ler(D-i) 

(D-l)(D-4) 



-17+2tt z 



AW + 



iorp-i) 8 

(D-l)( J D-4) + 3 



a 2 ^ + 0(D-4)|+0(e 4 ).(146) 



Note again that we have dropped terms which vanish for D = 4, and also 
terms which fall off at late times, relative to the overall factor of a 2 . The 
only infrared logarithm in (11461) is precisely the one predicted at order e 2 by 
the stochastic analysis [33] . 



5 (n\<p*<p\n) 

It might seem that computing this VEV requires an analysis as extensive as 
what we have just done. However, it is possible to get this result from the 
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previous one by making use of the exact Heisenberg equation of motion for 
the scalar field operator, 



-9 v ' 2 



(147) 



Next act the scalar d'Alembertian on the gauge invariant product of tp*(x) 
and ip(x), and substitute fll47[) . 

1 



-g 



+2g^(D^)*D u <p + <p*-—=D lx (^g-g^D 1/ tp) , (148) 

= 2g^(D^yD uV + + t^-(<p V) 2 • (149) 

l + oZ 2 l + oZ 2 

Because it is valid to use the equations of motion inside functional integrals 
of gauge invariant operators [51] . we conclude, 



l+SZ 2 



+ 



n 



Q) + -^—(Q 
1+SZ 2 



(<p*<p) 2 n) .(150) 



Finally, recall that 5£ ~ e 2 , 8Z 2 ~ e 2 and 5 A ~ e 4 . Hence the result we 
require at one and two loop orders is, 



Q) = 2g" v (Q (Pu<p)*D v ip n) + 26£R(n <p*(p fi) + 0(e 4 ) . (151) 



To the order we are working, the last term in (j!5ip is just 2R = 2(D — 
1)DH 2 times the product of the one loop conformal counterterm (|65|) with 
the coincident scalar propagator, 

28£R(n\(p*(p\tt) = (-e 2 2 J D 7 (0) + 2Ae 2 H 2 5^ n ) (A(0) + 2Hn(a)) + 0(e 4 ) . 

(152) 

The first term on the right hand side of (11511) is just twice the trace of our 
result (11461) from the previous section, 



2<T (ft {Dn<p)*D v <p 



n) = -2D 



1-5Z 2 



D-l 
D 



kH 2 



+ e 2 2D l{ 0){A(0) + 2kHa)} + ^^r(^) {^ + 80} 



(An) 



2 > ^D-A 



2 tt2D — 4 



e 2 H 



losrp-i) 

{Att) d \{D-1){D-A) 



120 + 16tt 2 +O( j D-4) +0(e 4 ) .(153) 
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Substituting these two relations in (11511) results in complete cancellation of 
the divergent infrared logarithms to the order we are working, 



<p cp 



n 



-2D 



1-SZ 2 



D-l 



D 

+2^H^ Q {A(0) + 2kHa)} + ^-^r(^) {^ + 80} 



kH 2 

e 2 J?% r 12D 



2 tt2D — 4 



+ 



e 2 H 



108Tp-l) 
{Att) d \{D-l){D-4) 



(An) 



•120 + 16tt 2 +O(D-4) \+0(e 4 ) .(154) 



In fact we can eliminate any infrared logarithms, at order e 2 , by choosing the 
finite part of the conformal counterterm to vanish, 



<p cp 



n 



= -2D[1-5Z 2 

<5&n=0 

e 2 H 2D-4j 108r(D-l) 



D-l 
D 



kH 2 



120 + 16tt 2 +O( j D-4) \+0{e 4 ) .(155) 



(4tt) d \(D-1)(D-A) 
Note that on a function of ln(a), the scalar d'Alembertian gives, 



□/(ln( C 



-H 2 lf"(Ha)) +(D-l)f (\n(a 



(156) 



If the d'Alembertian of such a function is a constant K then we can recon- 
struct the function up to an integration constant, 



□/(ln(a)) = K 



/(ln(a)) 



K\n(a) 



(D-l)H 2 



+ constant . (157) 



We see from this and (11551) that choosing <5£ fin = results in the VEV of 
<p*(x)<p(x) possessing no In 2 (a) contribution at order e 2 . It is interesting to 
note that 5^ n = is also the unique choice which results in there being 
no significant late time corrections to the scalar mode functions at one loop 
order [52] . 

At order e° the VEV of <p*(x)<p(x) is just the coincident propagator, and it 
contains a single infrared logarithm. Hence the leading logarithm correction 
for this VEV at order e 2 would contain two infrared logarithms. We have 
just seen that choosing <5£ fin = causes this leading logarithm correction to 
vanish. That is another key prediction of the stochastic formalism 
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6 Discussion 



We have used dimensional regularization to compute the one and two loop 
VEV's of two gauge invariant operators in SQED. Our results fl 146j) and (11551) 
confirm two predictions of the stochastic analysis [33J: 

• That the leading log result for the coincident kinetic term is, 



n 



LL 



H% v \ 3,3 e 2 



16vr 2 2 D-A Att- 



ln(a) + 0(e 4 ln 2 (a)) 1. (158) 



That setting the conformal counterterm to, 




' + J + 0(D-A)\, (159) 



D-A 2 

(which corresponds to S^ n = 0) results in the coincident scalar norm 
having no leading logs at order e 2 , 



n 



^(xM X )\n) hh = + x £- ln(a) + 0(e 4 ln 2 ( 



Air 2 An 2 



(160) 



Two additional points deserve comment concerning the leading logarithm 
approximation. First, we saw in Yukawa theory |32j that leading logarithm 
corrections to the VEV's of passive fields can harbor ultraviolet divergences. 
From (11581) we see that the same can be true for leading logarithm corrections 
to the VEV's of differentiated active fields. In both cases the reason is that 
the ultraviolet cannot be ignored in any field which fails to contribute an 
infrared logarithm. By contrast, the VEV of undifferentiated active fields 
such as (11601) must be finite at leading logarithm order. Note that this would 
be true no matter what choice had been made for <5£fi n . 

Our second observation is that choosing (11591) prevents significant late 
time corrections to the scalar mode functions |52j at one loop. This seems to 
be the unique renormalization prescription which suppresses quantum secular 
effects as much as possible at order e 2 . However, there is no way to prevent 
significant corrections to the photon mode functions at one loop order [H |£l 

nig. 
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Although the principal application of this exercise has been to provide 
"data" for checking testing a leading- log resummation of SQED [33], the cal- 
culation is not without interest in its own right. Because both the operators 
whose VEV's we computed are gauge invariant, it would not have mattered 
which gauge we used. However, we found that working in Lorentz gauge (1291) 
greatly simplified the one loop self-mass-squared. It also permitted the two 
loop scalar kinetic operator f)158p to be expressed in terms of a single-vertex 
integral, rather than the 2-vertex integration that seems to be required for 
the leftmost diagram of Fig. 3. 

Another significant technical advance (for which see the Appendix) is 
that we have worked out procedures for integrating functions of y(x; x') by 
extracting covariant d'Alembertians. This has advantages over the technique 
of breaking y(x; x') = H 2 aa'Ax 2 up into factors and then extracting powers 
of the flat space d'Alembertian d 2 = rj^d^d^. The older technique was em- 
ployed in all previous computations of this sort [51 El El [121 EE3 ESI E3] ■ 
Although it produces correct results, one must work through many tedious 
cancellations between spurious infrared logarithms from the D-dependent 
powers of a which reside on ultraviolet divergences, and equally spurious 
infrared logarithms from the nonlocal, finite terms. The new technique orga- 
nizes the calculation so that these spurious infrared logarithms never appear 
in the first place. 



7 Appendix: Tables SI HI and [5 



The purpose of this appendix is to describe how to evaluate certain integrals 
of the form, 

j d D x'a' D {f{y ++ )-f{y + _)} and j d D x'a' D - 1 {f{y ++ )-f{y + _)} . 

(161) 

Section 3 has already discussed our method for segregating ultraviolet diver- 
gences on to delta functions. Briefly, the procedure is to extract covariant 
d'Alembertians using the identity, 



(a-l)(a-±)H 2 ^y 



□ /4\a-l fOL— D\ /4\a-l 

■(«) + (; 



D 



V«^y. (162) 



One continues applying this until the singularities l/y a 1 are integrable for 
D = A. At this stage the denominator (a— y) vanishes for D = 4. To segregate 
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this divergence on a local term one adds zero in the form, 



(4n)-H- D i5 D (x-x') □ , 4 >.#-i D ,D w 4 xf-i 



One then takes D^4 in the nonlocal terms. 

The results for the two potentially divergent integrands we require are, 

/ 4 \D-2 _ 2 (47t)TH~ D i5 D (x-x') □ f/ 4 x ' 

W ~ (i>-3)(d-4) r(f-i) 7=T~ " ^IV n( ^~ , 

+ 2 (^)Kx)-^ +0(D - 4) ' (165) 

/ 4 2 □ / 4 \£>-2 2 / 4 \D-2 

The +— results follow from these by dropping the delta functions and re- 
placing y ++ everywhere with y + _. 

Of course it is quite simple to evaluate integrals of delta functions! The 
finite, nonlocal terms require an additional partial integration to remove 
factors of l/y. The two identities we need are, 

4 

y 



= #{|^(!)-Kf)}+3K|)-2- (168) 



These d'Alembertians are extracted from the integration over rr' M and then 
acted after performing the integrals. Because the integrals can only depend 
upon a, the following identities are useful for acting the d'Alembertians, 



°(M.,).^, °(540)-2J^_I, (170) 

□ / 2 . .x fil . . o □ /ln 2 (a)x 21n 2 (a) 21n(a) 2 



/J 2 



(ln 2 (a)) = -61n(a)-2 , " (f^) = f^_£i^_£ . (m) 
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We will shortly describe how to perform the finite integrals which make up 
Tables S] and E] but it is best to first complete the discussion of the potentially 
divergent integrals of Table [31 Combining (11651) with (I167P and (I168p . we first 
express the integral as a divergent part plus a sum of differentiated finite 
integrals. These finite integrals are performed using Tables H] and O then 
the derivatives are acted using relations (ll69H17ip . The last step involves 
expanding some terms on the divergent part to obtain a simpler answer for 
tabulation. The steps for the bottom left entry in Table [3] are, 



d D x'a >D - 



2/4 



D-2 



z(47r) 



d A x'a' 4 



2 V 4 



(D-3)(D-4) H D T(^-1 
-(+-)} + §~ 2 



6 ln(^±) -(+-)}+/ dVa' 4 { 6 ln(^±) 



i(47T 



• + 
2 



J aVa /4 |ln 2 (^±) 
>+0(D-4), (172) 



# D r(f 



■1){(D-3)(D-A) 



+ 



H 2 

i(47T 



±ln 2 (a)-fln(a) + f-^ 
li 



9 



- ln(a) + 
2 



+ 



^ln 2 (a)-|ln(a) + |-^ 



18 



11 
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ln(a) 

■ ' 6 



§ln(a) + 

' ' O(D-A) 



H D T(§-l){(D-3)(D-4) 



+ 2 + 0{D-4) 



i(47T 



H D T(%-1) [D-A 



O(D-A) 



(173) 
(174) 
(175) 



The analogous a term gives 



4 \D-2 



i(4% 



— -l 

2 a 



J (D-3)(D-4)H D T(±-1) 
W J d ^i^)M"f) - (+")} + £ /^{ln 2 (^ 
-61n(^) -(+-)} + / rfVa' 3 {61n(^±)-(+-)|+0( J D-4) , (176) 



Z'(47T 



2a- 1 



J ff Z) r(f-l)\( J D-3)( J D-4) 



# 4 



JM«)-S 

31 | 7T 2 1 



In 2 (a) , ln(a) 
4a a 



+ — - - + 

24a ~ 6a 6 ' 



ln(a) 
2a 



4a 
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+ 



H 2 



|ln(a)-^ 



In 2 (a) , 21n(a) 
2a ' a 

, 31n(a) 3_ 

^ a 2a 



+ 



31 | 7T 2 

12a 3a 



31n(a) 3 
1 — H 

a 2a 



i(47T 



2a" 



iW(f-l)\ (£-3)00-4) 



+ 0(Z>-4) 



i(47T 



2a 



-i 



i/ D r(f-i) [d-4 



0(D-A)\. 



0(D-4)|,(177) 
(178) 
(179) 



The top two entries in Table [3] follow from using relation (I166p on (11751) and 
(11791) . respectively. 

We turn now to the finite integrals in Tables H] and [5j First note that 
relations (I167P and (11681) allow us to express the results for f±(x) = 1/x and 
/ 2 (x) = ln(x)/x in terms of / 7 (x) =ln(x) and fs(x) =ln 2 (x), 



-I 



180) 



+3 { /7 (^)- /7 (^ 



181) 



So one uses these relations, and the d'Alembertian identities (1169111711) . to 
derive the first two entries on Tables H] and [5] from the bottom two entries. 
The procedure for integrating fr(x) and fs{x) is straightforward: 

• Combine the ++ and H — terms to extract a factor of i and to make 
causality manifest, 



>n(^)-ln( 



y+- 



27ri0 (A77 - 



ln2 (x)" ln2 (^f) = 47 ^K 

Here we define A?? = r] — rf. 

Make the change of variables r = x! 
integrations, 







x — 


^ ) 






\x— 


x'\) 



(182) 
(183) 



x and perform the angular 
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J d 3 x'9(Ar]-\x-x'\^F(-y(x;x') 

= A-Kd(Ar]) J drr 2 F(aa'H 2 (Ari 



2 2 

r 



(184) 



Make the change of variables r = Ar] ■ z and perform the integration 
over z. 

Make the change of variables a' = —l/Hrf and perform the integration 
over a'. 

Discard terms which fall off at late times with respect to 1 for Table HJ 
or to 1/a for Table [5j 



For the case of fs(x) integrated against a' 4 , the steps are, 

= 4m x 4nf\'a' 4 Ar] 3 f 1 dzz 2 \\ii( y aa) + 2\ji( y HAr])+\n 
hi Jo I 



1-z' 



16n 2 i r a 



16n 2 i f 1 



— In 

3 



^^ 2 f ^ 8 , / x 7 ^ 

, i - m a In (a H 

H A 6 K J 9 K 'A 9 



1 1 



(185) 
(186) 



a 2a 2 9 



13 



1 00 w _— n— 1 

^ L: 



'a° 



,-n-2 



ln(a) + — 

V ' 9a 

,. — n— 3 



12a 2 27a 3 ~^3n 2 n(n+l) n{n + 2) 3n(n+3) 
16vr 2 z J 1 
H 4 I 6 



,(187) 
(188) 



Note that it is only at the last step for which it matters whether the measure 
factor is a' 4 or a' 3 . 

The procedure is almost the same for fs(x) =ln(x)/(l — x) 2 and fi(x) = 
hi(x)/(l—x), with two exceptions. First, one extends the range of z and 
makes the change of variable, 



x = 



1 


r 1 


1 




-->} 


2 


■af~ 


a 





2X-4- 1 

q' q 

J__ 1 

q' q 



(189) 
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This carries the factor of At] 3 times the z integral to, 



(i-i) 3 \\zz*f{ 
v o ' a 1 Jo \ 



>i 1 1 
1 + aa 



2A-Z 1 



dx(2x 



/ / 1 1 
aa x[ — I — 



.(190) 



The second exception is that one changes variables in the temporal integra- 
tion from a' to a = a'/a. For fs(x) integrated against a' 3 , the steps are, 



d 4 x'a' 3 < 



' V++ 

^4 



2iti x 47r 
_ H 1 



da a / dzz" 



1 



[i+^'6-z) 2 (¥) 

^ [2x-(4 + i )l 2 



(191) 
(192) 
(193) 
(194) 
(195) 
(196) 



Note that the a' 4 result (on Table Hj) would follow from simply multiplying 
the integrand of (j!94p by aa. 

The most difficult reductions are those for fs(x) = ln 2 (x)— 2 ln(l— x) ln(x) 
and fe(x) = x ln 2 (x) — 2x ln(l — x) ln(x) — 2 ln(x). These combinations of fac- 
tors were chosen because they arise in performing the various integrations 
and multiplications, and because cancellations between the different factors 
prevent significant contributions from the limit at a = 1. We carry out the 
same reduction as for fz(x) and f±(x) on each factor separately, up to the a 
integration, and then combine them to take advantage of the cancellations. 
The two factors in /s(x) times a' 4 produce, 




In 2 fx) 



) =► - 1Fr j 1 — (l-«) 3 {--ln(«) + -ln(l-a)--},(197) 
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-21n(l-x) \n(x) 



16n 2 i [ x da < 1 



h 4 



9 9 



(198) 



These two factors sum to give the full result for f$(x) times a' 4 , 



~1F 



!l da { 

a \ 



2 + 



a 



ln(a) + _(l- a ) 3 hi(l-a) + --^ 



16n 2 i \ 1 vr 2 
H 4 I 6 ~ ~9~ 



+ 0(^)1. (199) 

Of course the a' 3 result follows from multiplying the a integrand by 1/aa. 

As might be expected, the reduction of fe(x) is the most complicated. 
The factor of xln 2 (x) times the a' 4 measure factor gives, 

, o, . 16n 2 i r 1 da , ,A 1 , , . 1 , . . 31 ] 
*b.»(*) ^ ^4 ^d-aJ^lnW-igMl-^ + ijo} ■ (200) 



The second factor produces, 

167r 2 i f 1 da 



-2x ln(l— x) In (a;) 



H 4 



r 1 da 
A a 2< 



1 a 2a 2 2a 3 a 4 a 5 



30 6 



6 30- 1 



31 11a a 2 a 3 11a 4 31a 5 
+ — - + - — + 



And the final term contributes, 

-2\n(x) = 

The sum of all three terms is 
16n 2 i r 1 



450 18 9 9 



16n 2 i r 1 da 



18 450 



H 4 



f L da a xS 
h * -3 (1 " a) ' 



ln(a) 
.(201) 

(202) 



H 4 



da< 



a a 2 a 3 
l_l 1 

3 3 15 



ln(a) 



1 c 1 9 9a a 2 1 

— ^(l-a) 5 ln(l-a)--^ + ^-^ + ^- , (203) 
lba z 15a 5 5 15 



16n 2 i f 5 7r 2 



H 4 24 18 



(204) 



As always, one obtains the result for the a' 3 measure factor from multiplying 
the a integrand by 1/aa. 
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